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ABSTRACT 



We present an update of the orbital ephemeris of the binary X-ray pulsar Her X- 1 and determine an improved value for the rate of 
orbital decay. In addition, we report the first measurement of the orbital eccentricity. We have analyzed pulse timing data of Her 
X-1 from X-ray observations by RXTE (Rossi X-Ray Timing Explorer) and INTEGRAL over the period 1996 - 2007. Accurate 
pulse arrival times were determined from solar system bary-centered photon arrival times by generating pulse profiles averaged over 
appropriately short integration times. Applying pulse phase connection techniques, it was possible to determine sufficiently accurate 
local ephemeris data for seven observation periods distributed over 12 years. Combining the new local T n /2 values with historical 
values from the literature we update the orbital ephemeris of Her X-1 to T nj2 = MJD 46359.871940(6) and P mb = 1.700167590(2) d 
and measure a continuous change of the orbital period of dP m b/dt = -(4.85 ± 0.13) x 10~"ss~'. For the first time, a value for the 
eccentricity of the orbit of Her X-1 is measured to be e = (4.2 ± 0.8) x 10~ 4 . 

Key words, stars: binaries: general - stars: neutron - X-rays: general - X-rays: binaries - X-rays: individuals: Her X-1 - Ephemerides 
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1. Introduction 

Hercules X-1 and Centaurus X-3 were the fir st two X-ray pul- 
sars d i scovered by the Uhuru satellite in 1972 dTananbaum et alj 
fl97l iGiacconi et ail fl973h . Since then they have remained 
among the most studied X-ray sources. Her X-1 is a close binary 
system consisting of an accreting magnetized n eutron star a nd 
the stellar companion HZ Her (first suggested by Liller 19721). a 
main sequence star of spectral type A/F (Crampton 197 4S). The 
mass of the optical companion is ~2. 2M„ dDeeteretaljl 19811) 
which places the system in between high and low mass X-ray 
binaries. The orbital period of the system is 1.700d. Due to the 
high inclination of the orbital plane (i ~ 85° - 88°) the X- 
ray source is regularly eclipsed by the optical companion for 
~5.5h once per orbit. The X-ray luminosity o f the source is 
Lx ~ 2 x 10 37 ergs~' for a distance of ~7 kpc (Reynol ds et alj 
fl997h .The spin period of the neutron star, represented by strong 
X-ray pulsations, is 1.237 s. 

The very first Uhuru observations of Her X-1 revealed the 
presence of a long-term ~35d periodicity which manifests it- 
self mainly through the alternation of so-called on (high X- 
ray flux) and off (low X-ray flux) states. The 35 d cycle con- 
tains two on states - the main-on (~7 orbital periods) and the 
short-on (~5 orbital cycles) - separated by ~4-5 orbital cy- 
cles. This periodicity in Her X-l/HZ Her is usually attributed 
to a counter-orbitally precessing tilted and warped accretion 
disk, the outer rim of which periodically blocks the line of 
sight to the X-ray emitting regions on the surface of the neu- 
tron s t ar (see, e.g., Gerend & Bovntonll976tlHow arth & Wilson 



119831: IShakura et alj 119991) . Roughly every ~ 5 yrs Her X- 
1 exhibits so-called anomalous low states where the X-ray 
flux is strongly reduced and the 35 d variability is seen only 
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marg i nally dParmar et alll985t IVrtilek etaDl994t iParmar et all 
ll999UStaubert et alj|2009l) ~ 

Here we report on the timing analysis of X-ray data ob- 
tained with RXTE (the Rossi X-ray Timing Explorer) and with 
INTEGRAL. Combining these results with historical data from 
the literature we refine the global orbital ephemeris of Her X-1 
and rate of the secula r decrease of its or bital period, which was 
originally reported bv lDeeter et al.l d!991l) . In addition, we were 
able to measure the eccentricity of the orbit for the first time 
(only upper limits had been reported up to date). 

2. Observations 

Her X-1 was repe atedly observed b y the Rossi X-ray Timing 
Explorer (RXTE) (Brad t et aDl993l) over a 10 years time frame 
(1996-2005). In addition, we make use of two extended ob- 
servations by INTEGRAL dWinkler et alJ 120031) in 2005 and 
2007. The observations provided high-quality X-ray spectral 
and timing data. For the timing analysis presented here we 
used da ta obtained with the RXTE Proportional Counter Array 
(PC A) dJahodaetal.ll 19961) in the energy ran ge 3-20 keV and 
with INTEGRAL ISGRI dUbertini et alj|2003h in the range 20- 
60 keV. Seven sets of data, all taken during a main on state, were 
long enough and sufficiently densely sampled such that a local 
ephemeris could be determined (the set of 2005 is from partly 
simultaneous observations by RXTE an d INTEGRAL, see also 
lKlochkovll2"007t iKlochkov et al.ll2008albl) . We note, that the ca- 
pabilities for relative and absolute timing required for the kind of 
analy sis presented here are well met by bo th RXTE dRots et al.l 
120041) and INTEGRAL dWalter et alj|2003l) . 

Details of the observations are given i n Table [p Fol l owing 
the 35 d cycle coun t ing m ethod adopted by Staub ert et al.l (l983) 
and I Staubert et al.l (120091) the seven cycles are numbered 257, 
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Table 1. Details of RXTE observations of Her X-l used for the timing analysis. 



^ ^ f.\if.\p* nil J 

jj u cycle nu. 


?S7 

Z.J / 








jaj 




j / j 


Source of observational data 


KA 1 r. 


KA 1 n 


KA 1 C. 


PYTC 


KA 1 r. 


KA 1 ryllN 1 £.LjK/\L 


11N 1 COKAL 


Date of observation 


July 96 


Sep 97 


Dec 00 


Dec 01 


Nov 02 


July 05 


Sep 07 


Start of observation [MJD] 


50290.027 


50704.912 


51896.319 


52243.060 


52600.061 


53574.715 


54346.088 


No. of RXTE/INTEGRAL orbits used 


18 


14 


33 


45 


29 


8/3 


1 


Elapsed time of observation [days] 


1.42 


2.43 


5.62 


3.52 


4.02 


8.16 


2.03 


Elapsed time of obs. [binary orbits] 


0.83 


1.43 


3.30 


2.07 


2.36 


4.80 


1.20 


Coverage of binary orbits [%] 


36.6 


10.0 


10.0 


30.7 


18.9 


37.0 


53.0 


Number pulse profiles 


139 


45 


61 


161 


133 


9/69 


27 


Mean integration time per profile [s] 


322 


276 


802 


580 


493 


420 / 2350 


3500 



1 The numbering of 35 d c ycles followslStaubert et al] dl983h : "Pulse profile counting" is applied, based on pulse profile shapes repeating with an 
average period of 34.98 d (Staubert et al. 20091). 



269, 303, 313, 323, 351 and 373. Here "pulse profile counting" 
is applied, as defined by Staubert et al. (2009) (meaning that the 
changing shape of the pulse profile is used to identify individual 
35 d cycles). The PCA data were mostly taken in "event mode", 
providing arrival times for individual events, and a few times in 
"binned mode", giving light curves (with time resolution, gener- 
ally 1/16 s, sufficient to do timing analysis for the 1.24 s pulsa- 
tion in Her X-l). 

3. Timing analysis of RXTE and INTEGRAL data 

3. 1 . Basic techniques to determine the pulse period 

We use two techniques for the determination of pulse peri- 
ods: e poch folding with x 2 search dHolt et al.lll976tlLeahy et al] 
1983b and pulse p hase connection dManchester & Tavlorlll977t 
Deeter et alj|1981l) . Both methods are generally well established 
and are appropriate in the case where the pulse period is al- 
ready roughly known. Our data are either "event lists" (the ar- 
rival times of single detector events) or "light curves" (contain- 
ing the numbers of events per second and the absolute times of 
the centers of those bins). 

With epoch folding, events are accumulated into bins of a 
phase histogram according to their phase with respect to an as- 
sumed trial period and a given start time (e.g., the time of the 
first event), thereby generating a pulse profile (e.g., Fig. Q] with 
our Her X-l data we generally used 128 bins). A search for the 
likely 'best' period is conducted by constructing a x 2 distribu- 
tion: for any pulse profile generated with an assumed trial period 
the value x 2 — £ («/ _ ( n )) 2 /cf is computed, where n,- are the 
number of events in bin i and (n) is the mean of all n,. x 2 is test- 
ing the deviation from a uniform distribution (which is expected 
for a trial period far from the true one). The x 2 distribution is 
found by plotting the x 2 values for a series of trial periods (se- 
lected to sample a sufficiently large period range with sufficient 
resolution). The 'best' period is assumed to be the one leading 
to the pulse profile with the maximum deviation from a flat dis- 
tribution. The x 2 distribution resembles a triangular distribution 
and can often be well fitted by a Gaussian, the mean of which 
may be considered the 'best' period (note that this is likely bet- 
ter than using the maximum of the x 2 distribution, especially 
in cases of limited statistics). The full width at half maximum 
(FWHM) of the x 2 distribution should be of the order FWHM 
~ P 2 /TOB, if P is the true pulse period and TOB is the total 
elapsed time of the observation. 

For pulse phase connection pulse profiles are produced, us- 
ing the best estimate for the pulse period (e.g., from ax 2 search), 
for integration intervals of appropriate length (leading to profiles 



with sufficient statistics). The time differences between profiles 
of adjacent integration intervals (measured either by making use 
of sharp features in the profile or by template fitting) are then 
identified with an integer number of pulse cycles. In this way, a 
list of unique cycle numbers and associated absolute times are 
generated, which can then be modeled. Necessary conditions for 
the applicability of this method are: (1) the shape of the pulse 
profile must be stable enough for the common details of differ- 
ent pulse profiles to be identified and (2) the time separations of 
the profiles must be small enough such that, with the given un- 
certainty of the assumed pulse period, any mis-counting can be 
ruled out. If these conditions are satisfied the pulse phase con- 
nection is by far the preferred technique, allowing very flexible 
modeling, especially if the pulse period is time variable, and with 
easy and accurate estimates of uncertainties. 

For a constant pulse (spin) period P s the expected arrival 
time of pulse number n is 



? n = t + n P s , 



(1) 



where to is a reference time at the start of the observation. For 
the case of non-zero first and second derivati ves of the pulse pe- 
riod, the arriv al times are given by (see, e.g.. lKellev et al.l ll980: 
NagasidUl) 



t D = to+nP s + ^n 2 P S P S + \ « 3 P 2 S P S + . 
2 6 



(2) 



Here f s is the initial pulse period valid at to- Modifications for 
pulse arrival time delays, e.g. because of binary motion of the 
pulsar, can easily be introduced (see below). Often Equ. <(2j 
is expresse d in the inverted form giving the pulse phase <j> at 
time ( t-to) dBlandford & Teukolsk^l976tlManchester & Tavlorl 



H977HKramer]l2008h 

1 ,1 , 

<p = 4>o + v (t - t ) + -v (t - to) 1 + -v(t-toY + ... (3) 

2 6 

with vq = l/P$ being the spin frequency at the reference time to- 



3.2. Analysis with previous ephemeris 

Our initial approach to the timing data of Her X- 1 obtained with 
RXTE was to determine accurate local pulse periods and to con- 
struct the corresponding pulse profiles. The scientific motivation 
was two-fold: (1) to complete our data base on the evolution of 
the pulse period with time which shows a long-term mean spin- 
up of about 10 ns per day. Superimposed to this trend, modu- 
lations are seen of relative spin-up and spin-down on different 
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Fig. 1. An example of PCA RXTE pulse profiles used to study 
pulse shifts due to the orbital motion of Her X-l. The profile is 
from 700 s integration time at MJD 52243.786 (35 d cycle no. 
313). The solid line shows the sharp decay after the soft trailing 
shoulder which was used as a phase reference (see text). 



time scales, some of which appear to be quasi-periodic and cor- 
related with the turn-on history o f the 35 d modulation due t o 
the precessing accretion disk, e.g.. IStaubert et al] (l2006i 12009). 
(2) With the spin periods found in this way, pulse profiles were 
constructed in order to study the systematic variations in pulse 
shape with phase of the 35 d modulation. This study has in fact 
led to support the idea of two different, but linked 35 d clocks in 
Her X- 1 : the preces sion of the accretion disk and free precession 
of the neutron star dStaubert et al.ll2Q09h . Further investigations 
with regard to the above two points are in progress. 

As a first step, all arrival times are translated to the iner- 
tial frame of the solar system barycenter. Then, the arrival times 
are corrected for the binary motion of the X-ray emitting neu- 
tron star around its optical companion, u sing the orb i tal ele - 
ments and ephemeris of Her X-l given by iDeeter et al.1 (Il99ll) : 
T„ /2 = MJED 43804.519980(14), P orb = 1.700167720(10) d, 
P orb = (-2.25 + 0.27) x 10- 8 dyr\ asinz = 13.1831(3)s, and 
eccentricity e = 0.0. Next, integration intervals are selected, gen- 
erally identical to the intervals of un-interupted observations (for 
RXTE the satellite orbits, for INTEGRAL the so called Science 
Windows), for which the local pulse period is determined by 
epoch folding/^ 2 search and the pulse profile is constructed. A 
typical profile is shown in Figure Q] Next, the phase connec- 
tion technique is applied to all profiles of the groups of observa- 
tions (Main-On states of Her X-l). Table[T]gives details for the 
7 Main-On states covered. This means that within each group 
every profile is associated with a unique pulse number and an 
absolute arrival time (in MJD). To find the latter we used a par- 
ticular feature in the pulse profile of Her X-l as a time marker, 
namely the "sharp edge" at the trailing edge of the right hand 
shoulder of the main pulse going into a flat bottom minimum. 
This flat bottom is generally well defined, extending over ~ 0.1 
in pulse phase over which the flux can be represented by a con- 
stant (for the example profile shown in Fig. 1 it is from phase 
~ 1.1 to ~ 1.2). A straight line is fitted to the "sharp edge" 
and the phase of the intersection between the straight line and 
the constant representing the minimum is determined. Since the 
"sharp edge" is only a few degrees from the vertical the system- 
atic uncertainty introduced by this procedure is quite small. The 
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Fig. 2. a) Delays of the pulse arrival time in Her X-l due to its 
orbital motion and theoretical sine curve for 35 d cycle no. 313 
(Dec 200 1). b) Residuals to a fit using the orbital parameters 
based on lDeeter et al.l (119911) . c) Residuals using the best-fit pa- 
rameters for our analysis of non-binary corrected data. 



overall accuracy has been found to be better than one bin for 
128 bins per pulse period, that is < 10 ms. Combining this with 
the known absolute time for phase 0.0, i.e. the reference time at 
which the folding of the profile was started, generally the time 
of the first event, gives the arrival time assigned to this profile. 
We have selected the "sharp edge" as the time marker, because 
it appears to be the most stable feature of the pulse profile of 
Her X-l, while most other features (including the main peak) 
change their shapes as a function of 35 d phase. We note, how- 
ever, that the timing results do not depend in any critical way 
on using the "sharp edge", consistent results have been obtained 
on the s ame data with template fitting using the complete pulse 
profile dKlochkovll2007t iKlochkov et aill2008bl) . The data sets 
used are short enough, such that the shape of the pulse profile 
does not change significantly. It is, of course, necessary that the 
accuracy of the assumed pulse period is sufficiently high, such 
that no counting errors occur when bridging the existing gaps 
between data sets. 

We like to mention here that the difference in energy range 
in which the data have been taken by RXTE (3-20 keV) and 
by INTEGRAL (20-60 keV) does not pose a problem for the 
current timing analysis, also not for the overlapping observations 
of RXTE and INTEGRAL of 2005 July. First, each data set was 
analyzed independently, leading to independent, but consistent 
solutions for the local T„p. values. Since it is known that the 
pulse profiles of Her X- 1 (and other pulsars) change with energy, 
one might expect that the absolute time of the "sharp edges" 
is energy dependent. Instead, we have found that the 2005 data 
sets of both satellites could be combined into a common set. 
Any shift between the two sub-sets is less than five milli-second 
(well within the general uncertainty), supporting the notion that 
the "sharp edge" is indeed a stable feature. 

Applying Equ. (0, the observed pulse number/arrival time 
data were fitted and the best fit values for to, P s and P s deter- 
mined (with P s = 0). A first set of T„/2 values was obtained 
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Fig. 3. Uncertainty contours for eccentricity versus omega from 
fits to the observations of 35 d cycle no. 313. The error contours 
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which were combined with historical T K /2 values from the liter- 
ature to establish a revised global solution from which updated 
values for P or b an d Arf> f° r an y given time could be calculated. 
However, we found that systematic sine-like residuals remained 
with an apparent period close to the orbital period, suggesting 
that the ephemeris which we had used to correct the arrival times 
for binary motion was not optimal. As an example, Fig. [2J3) 
shows the observed residuals for 35 d cycle no. 313 (Dec 2001). 
We therefore decided to perform a timing analysis (starting with 
barycenter-corrected, but not binary-corrected arrival times) and 
to determine the local ephemeris for all of our data sets. 

3.3. Analysis of non-binary corrected data 

In the same way as described in Section l3~2l the bary center (not 
binary) corrected arrival times are used to generate sets of pulse 
profiles and the corresponding lists containing pairs of pulse 
numbers and absolute arrival times. These lists were then taken 
for fits with equ. (f2]i which was extended by the additional term 

+ a sin i cos[2^ (f - T^/Porb] (4) 

in order to describe the arrival time delays due to the chang- 
ing position of the X-ray source along the binary orbit (as- 
sumed to be circular), a sin / is the projected radius of orbit in 
light-seconds and T n /2 (= Tgo) in MJD is the time at which the 
mean orbital longitude of the neutron star is 90°, correspond- 
ing to the maximum delay in pulse arrival time, or zero pulse 
period derivative (in Her X-l essentially ide ntical to the center 
of ecl ipse). For these fits a sin i = 13.1831s dDeeteret al.lfl98ll 
119911) and orbital periods P or b (determined for each data set from 
the global solution found earlier) were kept fixed. The fit param- 
eters, determined for all data sets, are 7^/2, to, P s , P s (P s = 0). 
For the data set of 35 d cycle no. 313 (Dec 2001), the improve- 
ment is evident in Fig . [2t) . 

3.4. Final fits solving Kepler's equation 

The results obtained so far allow us to establish a global solu- 
tion with quite accurate values for P or b and P or b for any time, 
including the times for all of our observations. Such a procedure 
is necessary since P 01 -b cannot be constrained from the individual 
data sets themselves. 



Keeping P or b fixed at values calculated from the above global 
solution, we have attempted to find a full orbital solution for 
each of our individual data sets by solving Kepler's equation 
dBlandford & Teukolskvl [1971 iNagasd fl989l) . in order to see 
whether we can constrain other parameters which were so far 
kept constant: the projected orbital radius a sin i, the eccentricity 
e, the longitude of periastron passage co (in case e is not zero) 
and the second derivative of the pulse period P s . 

The first parameter tested was a sin i. Our best data set (cycle 
no. 313, Dec 2001) gives a sin / = 13.1831(5) s. For all other data 
sets the values are consistent with this (within uncertainties). A 
weighted average of all values lea ds to a sin / = 13.18 3 1(4) s . 
The same value was found by iDeeter et al.l d 198 ll 1 1 99 lb : 
a sin/ = 13.1831(3) s (also a weighted average of several sets of 
observations). So we have used this value of a sin i as a constant 
input parameter. We note here, that we had wondered about our 
inability to reduce the uncertainty of this value, despite the fact 
that the number of photons we have collec ted from Her X-l i s 
about 20 times larger than that collected bv lDeeter et al] d 198 lb - 
We suspect that this may be due to the intrinsic variability of 
the pulse profile which introduces additional scattering in pulse 
arrival times, limiting the achievable accuracy independent of 
photon statistics. 

The second parameter tested was the eccentricity e. To our 
surprise, the fits of three data sets (35 d cycle nos. 257, 313 and 
323) yielded formal values for e larger than zero. For our best 
data set (cycle no. 313) a statistically significant value is found: 
e = (4.2 ± 0.8) x 10~ 4 , the corresponding value for the longitude 
of periastron is a> — (96 ±10) deg Cf^ in + 1 -0, innermost contour 
of Fig. [3]). This result suggests a significance of little more than 
5 sigma. In order to furth er investigate the statistical situation, 
we have applied an F-test (Bevington 1969): the reduction in^ 2 
when adding e and to as additional free parameters is 27.2 at the 
final level ofx\ = 156.6 with 156 degrees of freedom (dof). This 
leads to F x = (Ay 2 /,Y 2 )/(A(dof)/dof) = 13.5, corresponding to a 
probability of < 10~ 4 that the improvement in^f 2 is simply due 
to statistical fluctuations. Further confidence in a significantly 
measured eccentricity is gained by combining data sets. We have 
combined the data sets of cycle nos. 313 and 323 into one com- 
mon set by correcting the arrival times for the local P s values 
and shifting the arrival times of cycle 323 to the absolute time 
frame of cycle 313, using our best information on P or b and P or b. 
The fit to the combined data set yields e = (3.8 ± 0.8) x 10~ 4 
and <d = (82 ±10) deg, consistent (within uncertainties) with 
the results from cycle 313 alone. Adding the data from cycle 
257 or cycle 304 leads again to consistent results, but no further 
improvement. Even though the combination of several data sets 
confirms and strengthens the results of the single data sets, we 
prefer to quote the results from cycle 313 alone as our measure- 
ment, since it results from the simplest and most straight forward 
analysis of our best data set. From here onward we have kept e 
and 10 fixed at the values measured for cycle 313. 

Finally, we attempted to measure values for the second 
derivative P\ of the pulse period P s . Only for one data set (cycle 
323) a statistically significant value is found, for all others P s is 
consistent with zero, for cycle 373 also P s is consistent with zero 
(see Table|2]i. 

All data sets have then finally been fitted, solving Kepler's 
equation and leaving the following parameters to be determined: 
^s> T„/2, to, P s and P s . The final fit parameters and their uncer- 
tainties are summarized in Table [2] The uncertainties quoted are 
1 <x (68%), determined from the projections of the (x 2 mia + 1-0) 
error contours onto the parameter axis of interest, allowing all 
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Table 2. Orbital parameters of Her X-l determined from the RXTE observations during 35 d cycles No. 257, 269, 303, 313, 
323, 351 and from INTEGRAL observations during cycles 351 and 373. 1 The uncertainties in parentheses are 1-sigma (68%) 
^min + 1-0) an d re f er to the last digit(s). The following orbital elements were kept constant: a sin i = 13.1831 s, omega = 96.0deg 
and e = (4.2 + 0.8) x 10~ 4 . The values for P or b are calculated using the best-fit orbital elements given in Table[3] 



35 d cycle no. 1 


P s M 


T n/2 [MJD (TDB)] 


/Vb [d] 


dPJdt [ss-'] 




dPJdt/dt [r 1 ] 


257 


1.237739612(25) 


50290.659199(20) 


1.700167399 


-(2.8 ± 1.0) 10" 


12 


0.0 


269 


1.237729884(35) 


50703.799869(50) 


1.700167379 


-(6.8 ± 2.5) 10" 


12 


0.0 


303 


1.237769831(12) 


51895.617167(25) 


1.700167321 


+(1.54 ±0.17) 10" 12 


0.0 


313 


1.237767113(6) 


52242.451312(10) 


1.700167304 


-(5.7 ± 1.5) 10" 


13 


0.0 


323 


1.237761690(4) 


52599.486440(12) 


1.700167287 


+(2.8 + 0.2) 10" 


12 


-(2.9 ±0.08) 10" 17 


351 


1.237758409(8) 


53573.682256(35) 


1.700167240 


-(7.0 + 0.8) 10" 


13 


0.0 


373 


1.237744750(60) 


54345.558195(80) 


1.700167202 


0.0 




0.0 



1 The numbering of 35 d cycles follows Staubert et all il983h : "Pulse profil e counting" is applie d, based on pulse profile shapes repeating with an 

average period of 34.98 d dStaubert et~ai1l2009h . 



other inter esting parameters to vary during the error contour 
evaluation (Lampto n et al.|[T976b . 

The values of T„/2 are significa ntly different from those cal- 
culated using the global solution of iDeeter et al.l dl99ll) . This is 
not surprising, giving the extrapolation over ~ 20 yrs and the up- 
dated value for P or b. In the next Section |4] we will determine the 
updated global solution for P or b and P 01 -b, using all available his- 
torical observations of the source combined with our new results 
provided in Table [2] 



4. Orbital ephemeris 

Using our new values o f T„/2 (Table |2]l and the historical data 
from previous missions (IDeeter et al.lll98lT 1199 it IWilson et al.l 
11994 and references therein) we determine an improved value 
for the time derivative of the orbital period P 01 -b and an up- 
dated orbital ephemeris of Her X-l valid for the whole ~36 yrs 
of observations since the discovery of the source. Modeling all 
T„/2 values (historical observations plus RXTE and INTEGRAL 
data) by a linear ephemeris (P or b = 0) results in a very poor fit: 
^red = 58.9 for 28 d.o.f. Thus, the linear ephemeris can clearly 
be rejected. Figure|4](upper panel) shows the residuals of T„/2 af- 
ter subtracting the linear part of the quadratic best-fit ephemeris. 
For the quadratic ephemeris (Port) ^ 0, shown by the solid curve) 
the fit gives = 1.8 for 27 d.o.f.. If the outlying dat a point 
around MJD -45000 (due to Tenma, IDeeter et all 1 19911) which 
deviates from the quadratic fit by more than 3cr is removed, . 
is reduced to 1.3. 

The best fit parameters for the quadratic ephemeris are listed 
in Table [3] P or b(0) is the orbital period at the reference time 
7*72(0). For any particular observation T„/2 and P or b can be 
found using the following formulae: 



1970 1975 1980 1985 1990 1995 2000 2005 



P«b(f) = Porb(0) + Porb(0) " (t - T„, 2 (0)), 



T nl2 {n) = T n/2 (0) + n P orb (0) + ^n 2 P olb (0) • P or b(0); (6) 

where n is the orbital cycle number (n — corresponds to 

IDeeter et alj d!991l) had also discussed an alternative model, 
representing constant orbital periods before and after a sudden 
change of the orbital period in 1983 (at MJD -45000). This ap- 
pears to be mainly driven by the outlying Tenma data point. We 
have performed linear fits to the data before and after 1983, they 
are shown as dashed lines in Figure f4)l. While the data set before 
1983 is marginally consistent with a linear function (^ rea = 1.8 




40000 42000 



44000 



46000 48000 
Time [MJD] 



50000 52000 54000 



Fig. 4. Upper panel: Residuals of T„/2 after subtracting the lin- 
ear part of the quadratic best-fit ephemeris. For the RXTE and 
INTEGRAL observations analyzed in this work 35 d cycle num- 
bers are indicated. The solid line represents the best fit quadratic 
ephemeris. The dashed lines correspond to fits by the linear 
ephemeris with a sudden jump of the orbital period around 1983. 
The lower panel shows the residuals to the best fit quadratic so- 
lution. 



for 18 d.o.f.), the data set after 1983 is clearly not: Xmd = 9.1 for 
8 d.o.f.). We, therefore, reject the idea of a sudden change of the 
orbital period in 1983. 



(5) 5. Discussion 



The value of P orb = -4.85(13) x 10~ n ss" 1 determined in Sect.il 
and its uncertainty are signi ficantly smaller tha n the corre- 
sponding values determined bv lDeeter et al.l (1 1 99 ll) (-6.16(74)x 
10 -11 ). In order to connect the derived value of the secular vari- 
ation of the orbital period with the physics of mass exchange in 
the system. IDeeter et all ( 1 1 99 lb had considered a very simplified 
model, in which the two stars comprising the binary are treated 
as point masses and only orbital angular momentum is consid- 
ered. In the case of conservative mass transfer one can obtain: 



+ 



M, 



opt 



3P 



(7) 
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Table 3. Orbital elements of Her X- 1 . The uncertain ties in paren- 
theses (68%) refer to the last digit(s). Note that iDeeter et all 
d!98ll) had measured a sin i = 1 3. 183 1 (3) lt-s (with a slightly 
smaller uncertainty). All other values are updated or new. 



7X0) [MJD (TDB)] 
/WO) [d] 
/WO) [d/d] 
a sin i [lt-s] 
eccentricity e 
omega [deg] 



46359.871940(6) 

1.700167590(2) 

(-4.85 ± 0.13) x lO" 11 

13.1831(4) 

(4.2 ± 0.8) x 10~ 4 

96.0 ± 10.0 



where M as and M opt are the masses of the neutron star and the 
optical companion, respectively, and M ns is the accretion rate 
of the neutron star. In this idealized case, using M ns = 1.4M , 
M opt = 2.2M , and the value of P orb = -4.85 x 1CT 11 s/s, 
we obtain a mass accretion rate of the neutron star of M ns ~ 
1.3 1(T 8 M Q yr _1 ~ 8 x 10 17 g s" 1 . Th is value is slightly smaller 
than that found bv lDeeter et al.l d!981b but still roughly an order 
of magnitude higher than the accretion rate necessary to provide 
the observed X-ray luminosity Lx ~ 2 x 10 37 erg s _1 (assuming 
that the radiative efficie ncy of accretion is ~10%). To explain 
this IDeeter et al.l (Il98ll) had suggested that the mass transfer is 
not conservative (M ns t -M opt ) and matter may leave the sys- 
tem in the form of magnetically channele d wind from the X-ray 
heated surface of the companion star (see lRuderman et ai] [l989) 
carrying away angular momentum. We argue, however, that the 
model used for conservative mass transfer is very approximate. 
It does not take into account the angular momentum of the ac- 
cretion disk and the accretion stream, their dynamic action on 
the stars as well as the shift of the center of mass of the sys- 
tem during the mass transfer, which are k nown to be significan t 
factors influencing the final result (see e.g. lD'Souza et al.l l2006). 
Furthermore, one might expect exchang e of angular momentum 
by spin-orbit coupling (Mars h et al.l2004l) . Thus, Equ. (01 is very 
approximate and may provide only a rough estimate (to the order 
of magnitude) of M ns . The discrepancy we observe, therefore, 
cannot rule out a conservative mass transfer scenario. 

The alternative to a continuous decrease of the orbital pe- 
riod, namely a sudden change of th e orbital period around 1983 
as discussed bv lDeeter et al.l ( T981I). had received observational 
support by the analysis of IStelzer et al.l d 1997b . who performed 
a fit of the timing data of Her X-l available at that time, finding 
that the model with a suddenjump provided a smaller^- 2 than the 
model assuming continuous decrease of P orD . Our current obser- 
vational result clearly rejects this conclusion (Sect. |4). This is 
reassuring in the light of several th e oretic al difficulties which 
had been discussed by Dee ter et al.l d 198 lb . Her e we add that 
from t he disk-stream coupling model described in lStaubert et al.l 
(2009) we expect a decrease of the mass transfer to the neutron 
sta r during Anomalou s Lows (AL), not an increase as assumed 
bv lDeeter et al.l dl98 lb . This is supported by the observed spin- 
down during most of ALSs dStaubert et al.ll2006b . Thus, we ar- 
gue that the decrease o f the orbital period of HerX-l/HZHer 
originally discovered bv lDeeter et al.l d 1 98 lb is due to a constant 
negative P 01D and is not confined to any short time interval (like 
a jump around 1983) and that the most probable explanation is 
a continuous mass transfer from HZ Her to the neutron star in a 
conservative or non-conservative scenario. 



6. Summary and conclusions 

The precise timing analysis of the RXTE an INTEGRAL ob- 
servations in combination with the historical data allowed us to 
improve the value for the secular chan ge of the orbi t al peri od in 
Her X- 1/HZ Her, originally reported bv lDeeter et ail dl981h . Our 
ab solute value of P nr b is somewhat smaller than that determined 
bv lDeeter et al.l d 198 lb . and the uncertainty is significantly less. 
We conclude that we cannot rule out the conservative mass trans- 
fer scenario. The data rule out the possibility of a suddenjump of 
the orbital period around 1983 as an alternative for a continuous 
decrease. The new value for P or b along with the values of T„/2 
and P 01 b (listed in Table [3} establish a new orbital ephemeris of 
Her X-l that can be used to calculate the orbital period, P orD , and 
the time T„/2 for any particular time. 

A further result is the first measurement of a non-zero eccen- 
tricity e = (4.2 + 0.8) x 10~ 4 of the Her X-l/HZ Her orbit. For 
most practical purposes, such a small deviation from zero may 
not be relevant. However, in searches for subtle effects, e.g. a 
precession of the pl ane of the orbit, possibly associated with pre - 
cession of HZ Her dDeeter & Bovntonll976l;[Deeter et aljl991h . 
such accurate values may be needed. 
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